Entanglement of electrons in interacting 

molecules 

Tina A.C. Maiolo* Fabio Delia Salaj Luigi Martinaf Giulio Soliani§ 



Abstract 

Quantum entanglement is a concept commonly used with reference to the 
existence of certain correlations in quantum systems that have no classical 
interpretation. It is a useful resource to enhance the mutual information 
of memory channels or to accelerate some quantum processes as, for ex- 
ample, the factorization in Shor's Algorithm. Moreover, entanglement is 
a physical observable directly measured by the von Neumann entropy of 
the system. We have used this concept in order to give a physical meaning 
to the electron correlation energy in systems of interacting electrons. The 
electronic correlation is not directly observable, since it is defined as the 
difference between the exact ground state energy of the many-electrons 
Schrodinger equation and the Hartree-Fock energy. We have calculated 
the correlation energy and compared with the entanglement, as functions 
of the nucleus-nucleus separation using, for the hydrogen molecule, the 
Configuration Interaction method. Then, in the same spirit, we have ana- 
lyzed a dimer of ethylene, which represents the simplest organic conjugate 
system, changing the relative orientation and distance of the molecules, 
in order to obtain the configuration corresponding to maximum entangle- 
ment. 
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1 Introduction 



Much attention has been recently devoted to the notion of quantum entangle- 
ment for its potential applications in quantum information theory P^. Indeed, 
quantum entanglement holds a fundamental role in quantum teleportation [5], 
superdense coding [3], quantum key distribution [J, quantum cryptography [5], 
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and has been proved to be necessary for an exponential speed-up of quantum 
computation over classical computation. 

For these reasons, a lot of efforts have been devoted to characterize entangle- 
ment and experimental measurements have demonstrated that it can affect the 
macroscopic properties of solids [jj. For example, Chen et al. analyzed in [7| 
entanglement of the ground states in XXZ and dimerized Heisenberg spin chains 
as well as in two-leg ladder by using the spin-spin concurrence and the entangle- 
ment entropy between a selected sublattice of spins and the rest of the system. 
He et al. showed in |S] the energy spectrum, pair-correlation functions and 
the degree of entanglement of two-electron states in self-assembled InAs/GaAs 
quantum dot molecules via all-bound-state configuration interaction calculation 
and compared these quantities in different levels of approximations. Buscemi 
et al performed in [5] the quantitative evaluation of the entanglement dynam- 
ics in scattering events between two indistinguishable electrons interacting via 
Coulomb potential in ID and 2D semiconductor nanostructures. They followed 
the idea introduced by Schliemann JU| and extended by other outhors ^J, [T^) . 
since it allows to use a density matrix in the spatial coordinates and to com- 
pute the von Neumann entropy. Such an approach is based on an analogous of 
the Schmidt decompostion for state vectors of two fermionic particles: through 
an unitary transformation the antisymmetric wave function is expressed into a 
basis of Slater determinants with a minimum number of non-vanishing terms. 
This number, known as Slater rank, is a criterion to identify whether a system is 
entangled or not, which involves the evaluation of the von Neumann entropy of 
one particle reduced density matrix. In order to exploit entanglement in quan- 
tum computation also more complex molecules are studied. For example, the 
researchers of Halm -Mcitncr Tnstitut are working on the realization of molecu- 
lar quantum computers whose quantum bits could be formed by fullerenes with 
a single nitrogen or phosphorus atom trapped inside. 

For this reason, it is important to investigate the properties of some molecules 
and the interaction among them. A study of this kind is carried out in this pa- 
per. Here, the starting point is the so-called Collin's conjecture ^S|> that is> the 
correlation energy in molecular systems is proportional to their entropy. This 
conjecture was confirmed by numerical evidence by Ramirez et al. an d taken 
up by Huang and Kais in |12| who interpreted the degree of entanglement as an 
evaluation of correlation energy. Entanglement is in fact a physical observable 
directly measured by the von Neumann entropy of the system. This concept 
is used in order to give a physical meaning to the electron correlation energy, 
which is not directly observable since it is defined as the difference between 
the total energy of a given molecular system, with respect to the one obtained 
with the Hartree-Fock approximation method. The Hartree-Fock method, in 
fact, is typically used to solve Schrodinger equation for a multi-electron atom 
or molecule described in the fixed-nuclei approximation (Born-Oppenhcimcr 
approximation) by the electronic molecular Hamiltonian and the calculation of 
the error due to this approximation is a major problem in many-body theory 
and a vast amount of work has been done on this subject |15|. 

In this work, we make a measurement of electron entanglement in two dif- 
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ferent examples of bipartite systems, as H2 molecule and the dimer of ethylene, 
where each hydrogen atom or each ethylene molecule, respectively, can be con- 
sidered a qubit. In this article, firstly, we decide the method we use to make 
a measurement of entanglement comparing the von Neumann entropy of the 
reduced density matrix S(pi) and the von Neumann entropy of the density ma- 
trix S(p) for the H2 molecule. After, fixed that the most convenient method 
to adopt is the second, (S(p)), we make a measurement of entanglement in a 
many-body system represented by a dimer of ethylene and we compare it with 
its electron correlation energy. Finally, using the Klein's inequality for a bipar- 
tite system, we define the interaction electron entanglement in order to study 
the degree of entanglement between the two molecules of the dimer. 

The article is organized as follows: in Sec. 2 we briefly resume what entan- 
glement is m3 (equivalently, when a state is separable or entangled and 
why we choose the von Neumann entropy to make an estimate of the degree of 
entanglement. In Sec. 3 we first analyze a dimer of ethylene considering it as 
a multielectrons system, calculate the correlation energy, and then compare it 
with the entanglement, as a function of the distance between the two molecules. 
In Sec. 4 we introduce a new quantity, the interaction electron entanglement, 
which is defined as S int — 5[2C 2 i?4] — 2S , [C 2 i?4], in order to obtain the entan- 
glement between the two molecules of the dimer. Then, we change the relative 
orientation and the distance between the molecules, in order to obtain the con- 
figuration corresponding to maximum entanglement. In this way, the system 
can be considered bipartite and each molecule can be seen as a qubit for an 
application to quantum computing. In Sec. 5 we explain some computational 
details describing the package we use and how we prepare the input for the 
compound systems that we analyzed. Finally, in Sec. 6 some comments and 
results are discussed. 

2 A measurement of entanglement 

In quantum mechanics one generally distinguishes between pure states and 
mixed states depending on the information one has on the state in which a 
physical system actually is [T8|. 

Pure states are usually introduced whenever one wants to describe a pysical 
system that is in a given state with certainty, and are represented by unit 
vectors \if>) of the Hilbert space associated with the system. Equivalently a pure 
state can be represented by a density operator p = which is such that 

Trp = Trp 2 = 1. 

Whenever one doesn't know exactly in which state the system is, but there 
is only a probability p, for the system to be in the pure state represented by 
then one tipically says that the system is in a mixed state and represents it 
by the density operator p = J2i Pil^i) = 12iPiPh where p, > 0, J2iPi = h 
and pi is the density operator associated with the pure state Contrary to 
a pure state, a mixed state is such that Trp = 1 ^ Trp 2 . 

The most important consequences of the mathematical structure of the state 
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space reveal in the study of compound systems. In fact, let us consider a physical 
system T made up by the two subsystems A and B whose states are represented 
by the elements of the Hilbert spaces J^a and J^b , respectively. Then, the sys- 
tem T is described by Mr = j^a ® -^b and a pure state of T can be represented 
by the unit vector S Mr- Let {\(f>i(A))} ie ^ and {\Xj(B))}jeN be orthonor- 
mal bases on the Hilbert space M'a and ^b, respectively. A pure state of 
,ff T = M'a ® -%?b can be expressed in the form \tjj) = J2i j c n \4>i(A)) <8> \xj{B)), 
where Y^i j \ c ij\ 2 = 1- Thus, if the state S J$?t can be written as a ten- 
sor product of an element of J4?a and an element of J4?b, as in the following, 
[tp) = \(j)(A))<g)\x{B)) = \(j)(A) , x(B)} , we say that is a product (or separable, 
or factorizable) pure state. A pure state of a pair of quantum systems is called 
entangled if it is unfactorizablc. According to a well known theorem due to von 
Neumann J^j, the vector \tfj) 6 J^t can be expressed in the so-called Schmidt 
decomposition: \tp) — 53i>i \fPi\ v i{A))\\{B)) . Then, \ip) is a product state if 
there is only one Schmidt coefficient different from zero (if z = 1 in the sum). 
On the contrary, if there are more than one non-zero Schmidt coefficients (if 
i > 1 in the sum) the pure state is entangled. 

As a simple bipartite system, let us consider the H2 molecule, whose Hilbert 
space can be represented by: 
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where ££ and represent the orbital and the spin space, respectively; in brack- 
ets we denote one of the two electrons while the index represents the dimension 
of the space. 

In the occupation number representation 20 [n\ t, ri\ J., n% ] \ n 2 I, ■ ■ ■ , n m t 
; n m I) the subscripts denote the spatial orbital index and with this notation 
let us introduce an orthonormal basis for each space c ^ m : 
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A pure two-electron state \^t) can be written in this representation as 

m m 

i*>=EX>».*4<4i°>> 0) 



a=l 6=1 



where |0) is the vacuum state, the coefficients u> a ,b satisfy, accordingly to Pauli 
exclusion principle the following requests: 
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and and c are the creation and annihilation operators of single-particle states, 
respectively, whose action on the vacuum state is 



c+c+|0) = \ab), a, b €{1,2,3, 4,... ,2m} 



3 ee |na T> 4 ee |ti 2 |) 



TO - 1 = | Tim T) m= \n m I) 

(5) 

The coefficients uj a ,b can be calculated by using the configuration interaction 
method |2J> In particular, for ff 2 molecule we use configuration interaction 
single-double (CISD) method that is limited to single and double excitations: 

|* CJSD ) = c ^o)+E«)+ E CfclO ( 6 ) 

ar a<6,r<s 

where |^o) is the ground state Hartree-Fock wave function, c r a is the coefficient 
for single excitation from orbital a to r, and c ' t is the double excitation from 
orbital a and 6 to r and s. In the occupation number representation, the CISD 
wave function is given by 

|*> = C0II1I2O3O4 . . .)+c?|0il 2 l304 ■ ■ • )+^|li0 2 3 l4 ■ • • )+c^|0i0 2 l 3 l4 ...) + ••• 

(7) 

where each "1" represents the presence of the electron and the subscript repre- 
sents the site where the electron is |2~3*] . 

In order to realize whether the state of the compound system is entangled or 
not, we have to construct the corresponding reduced density matrix. Indeed, the 
vector \ip) represents a product state of the compound system if and only if the 
corresponding reduced density operators represent pure states (equivalently, \ip) 
is an entangled state if and only if the corresponding reduced density operators 
formally represent mixed states). Hence, if we want to measure the degree of 
entanglement we can calculate the degree of mixedness of the corresponding 
reduced density operator (or matrix) [21] . This is a consequence of the fact that 
an entangled state gives rise to correlations whose information is lost if we only 
observe the subsystems, as it results from the fact that we effectively have a 
mixed state. 

Of particular significance for describing randomness of a state and mixedness 
of density operators is the von Neumann entropy, defined as 

S(p) = -Tr(plog 2P ). (8) 

It is important to note that p is a non-negative trace class operator, while S is 
not. In analogy with classical entropy, 5* measures the amount of randomness in 
the state p. More precisely, the entropy S is zero if and only if the state is pure, 
and it is maximized when S{p) — logid, where d is the dimension of the Hilbert 
space. In other words, the more mixed is the reduced density operator, the 
more entangled is the original state and this result can be seen as a justification 
for the use of entropy as a measurement of quantum entanglement. Thus, the 
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entanglement of a pure state of a pair of quantum systems can be obtained 
as entropy of either member of the pair, consequently, in order to calculate 
the degree of electron entanglement in the H 2 molecule we have to calculate 
the von Neumann entropy of the density matrix reduced with respect to an 
electron. Starting from the density matrix p = |v]/ CISD )(\I> CISD |, we obtain the 
reduced density matrix by making the partial trace relative to all the occupation 
numbers except n\ ] (see Appendix A) 

n niJ,nj I 1 



{ 1 



1 < 3 < : 



= E ( n i l, n 2 ^n 2 I,..., n rn ^,n m I \p\nx I, n 2 ^,n 2 I,..., n m ^n rn I) 



rii T= 0,1 
nj J = 0, 1 



/ u2*+l|2 , y^™- 1 I 2i+l,2j+2i2 

n CISD _ I l^i=\ l c l I + l^i,j=l l C l,2 I 

" 1 I n I 1 2 i v^ m_1 l„2i+2|2 



o \co\ 2 + ET=V 
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(9) 

The von Neumann entropy of the reduced density matrix p CISD represents the 
degree of entanglement: 

S(P? ISD ) = -Tr(p? ISD log 2 p? ISD ) = (10) 

m—1 m—1 m—1 m — 1 

= -(E i c " +1 i 2 + E | C ?: 2 +1 ' 2i+2 i 2 )^ 2 ( E i c " +1 i 2 + E i^2 +1 - 2j+2 i 2 )+ 

z z— 1 z i— 1 

m—1 m — 1 

-(ki 2 + e i^ +2 i 2 )^ 2 (m 2 + e i c " + T), (ii) 
i=i i=i 

In order to make a generalization of this method to multipartite systems, we 
decide to measure the degree of entanglement as the von Neumann entropy of 
the density matrix corresponding to the compound system: 

= -Tr(^plog 2 p^j = ~(^2ailog 2 ai + ^[3 l log 2 [3^, (12) 

where on and /3, are the so called Natural Spin Orbitals (the eigenvalues of 
p). We compare the behaviours of the entanglement obtained with the CISD 
coefficients method or with the Natural Spin Orbital method, calculated as a 
function of internuclcar distance R. 
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Figure 1: The figure shows 
the degree of entanglement 
in the Hydrogen molecule 
calculated as the von Neu- 
mann entropy of the reduced 
density matrix with respect 
to an electron (S(pf ISD )) 
and as the von Neumann en- 
tropy of the density matrix 
(S(P)). 



R(A) 



In Figure 1 we show that the entanglement calculated as the von Neumann 
entropy of the reduced density matrix S{p c ( ISD ) has a similar behaviour as the 
entanglement calculated as the von Neumann entropy of the density matrix S(p) 
of the compound system. In particular, at small intcrnuclear distances, the two 
curves representing the entropies overlap, while outside the bond region, they 
give different results though their behaviours are similar. Since we are interested 
to a qualitative course of entanglement and to the search of its maximum value, 
we neglect the above differences and choose to adopt the Natural Spin Orbital 
method in order to investigate more complicated systems. 



3 Entanglement as a measurement of electron 
correlation 

Our main aim in this Section is to show that entanglement can be interpreted 
as a measure of the electron correlation E3- Entanglement, in fact, 

is a physical observable directly measurable by the von Neumann entropy of 
the system. On the contrary, the electron correlation energy is not directly 
observable, since it is defined as the difference between the exact total energy 
of a given molecular system, with respect to the one prescribed by the Hartree- 
Fock approximation method, that is 

E c = E^ ct - E HF . (13) 

The correlation energy is the energy recovered by fully allowing the electrons 
to avoid each other and Hartree-Fock method improperly treats interelectron 
repulsions in an averaged way. In other words, the exact wave function for a 
system of many interacting electrons is never a single determinant or a simple 
combination of a few determinants. The corresponding error is due to the 
correlations that are the analogue of the quantum entanglement in separated 
systems and are essential for quantum information processing in nonseparated 
systems. 

Since in ^5] t ne authors discussed the formation of H2 molecule, showing a 
qualitative agreement between the entanglement and the correlation energy as 
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r ( a ) Rotation angle (6 Q ) 



Figure 2: A comparison between entanglement and correlation energy in the 
dimer of ethylene, calculated by changing the mutual distance of the molecules 
(left panel) with molecules in a face to face orientation and by changing the 
relative orientation (right panel, the molecules are parallel at R=3.5 A). 



functions of nucleus- nucleus distance, let us generalize the argument proposed 
above to more complex molecules, as a dimer of ethylene. Set the distance 
between the molecular planes, we calculate the electron correlation energy of 
the dimer by the procedure described in Sec. 5. For the same configuration, 
we calculate the degree of entanglement with the Natural Spin Orbitals method. 
We repeate the same procedure changing the distance between the planes or, 
once we have fixed this distance, we make a rotation as in picture and the results 
are in the first and in the second panel of Figure 2 respectively. 

As we can see from the above picture, the entropy (S) and the correlation 
energy (E c ) have a comparable behaviour, although we use two different scales 
to represent these two quantities. In fact, when we plot them as functions of 
the plane distance, the two curves quickly decrease till R = 4.5 A, moreover for 
bigger distances, both of them decrease more slowly. When we set the distance 
at R = 3.5 A and change the mutual orientation of the molecules on their planes, 
we can see that both the entanglement and correlation energy decrease when 
the rotation angle grows up. The reason of this fact is that when a rotation of a 
molecule on its plane is made, the bigger is the rotation angle, the small becomes 
the superposition of the external orbitals of the two molecules. Consequently, 
the electron entanglement and the correlation energy decrease as it is shown in 
Figure 2. 

The above reasonings suggest us to interpret the entanglement as an effi- 
cient instrument to evaluate the electronic correlation energy, not only for H2 
molecule ^21; but also for some bigger systems. 

4 Results: interacting molecules 

The study of the degree of electronic entanglement that we have made in this 
paper is a useful resource for quantum computers whose input states are con- 
structed in order to be maximally entangled. In this way, it is possible to obtain 
an exponential speed-up of quantum computation over classical computation. 
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Therefore, the main aim of this work, is to analyze a dimer of ethylene, which 
represents the simplest organic conjugate system, in order to find the config- 
uration corresponding to a maximum degree of entanglement between the two 
molecules. Thus we have to consider the compound system as a bipartite sys- 
tem where each molecule can be seen as a qubit and, consequently, calculate the 
degree of the entanglement due to the interaction of the molecules only, thus 
neglecting the internal interaction of each molecule. 

In order to realize our goal we utilize the well known fact in quantum chem- 
istry that the correlation energy between the two molecules of ethylene can be 
obtained as the difference between the correlation energy due to the interaction 
of all the electrons in the compound systems and the correlation energy of the 
electrons in each molecule, i.e., 

E l L nt = E C [2C 2 H 4 ] - 2E C [C 2 H 4 ]. (14) 

Let us now observe that when the distance between the two C 2 H 4 molecules is 
infinite, we can consider them as two separated subsystems of the compound 
system dimer, which we denote by 2C' 2 H 4 , since they are uncorrelated. In this 
case, we have that p[2C 2 H 4 ] = p[C 2 H 4 ] ® p[C 2 H 4 ], thus, by using the definition 
of von Neumann entropy, we have that S(p[2C 2 H 4 ]) = 2S{p[C 2 H 4 \). In the 
general case, i. e., for finite disctances, the Klein's inequality holds, that is, 

S(p[2C 2 H 4 ]) > S( P [C 2 H 4 }) + S( P [C 2 H 4 }) = 2S(p[C 2 H 4 ]). (15) 

We thus define the interaction electron entanglement (<%„*)) as 

S (int) = S(p[2C 2 H 4 ]) - 2S(p[C 2 H 4 ]). (16) 

Then, we make a study of the degree of the interaction electron entanglement 
in the dimer of ethylene by changing the relative orientation and distance of the 
molecules. The results obtained are reported in Figure 3, where it is shown the 
degree of interaction electron entanglement as a function of the rotation angle 
(in (a) and (b)) and as a function of molecular distance (in (c)) calculated in a 
face to face configuration. We can see that, for R < 3.1 A, Sj nt grows up with 
the increasing of the molecular distance and the maximum value of this function 
moves from a rotation angle = 50° to 6 — 30° and in general, the closer are 
the molecules the smaller is the rotation we would make to obtain the maximum 
value of Sint- Instead, for R > 3.1 A, Sint decreases with the increasing of R, 
and becomes less sensible to the rotation angle since the molecules are more 
and more uncorrelated. In order to confirm this behaviour, we plot Si nt as 
a function of R for different torsion angles, in the Figure 3c. It is evident 
that the configuration obtained setting R = 3.1 A is the critical configuration. 
Indeed, for smaller distances, the entropy increases, while for bigger distances, 
it decreases with the increasing of R. Moreover, for R = 5.5 A, the correlation 
between the two molecules is approximativcly zero indipendently of the rotation 
angle. 
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Figure 3: Interaction electron en- 
tanglement for the dimer of ethylene 
as function of rotation angle, calcu- 
lated for different molecular planes 
distances (R) in a face to face con- 
figuration (in (a) and (b)). The 
same results are plotted as function 
of molecular distances for different 
rotation angles (in (c)). 



5 Computational details 

In order to calculate the CISD coefficients of the wave function, or equivalently, 
the entries of the reduced density matrix p^ ISD and to found the eigenvalues 
of the density matrix, the so called Natural Spin Orbitals, we use the package 
Gaussian 03 |2B]- In the calculation concerning the dissociation of the Hi 
molecule, we expand the wave function with a Configuration Interaction Single 
Double method. For this system, all the possible excitations are the single 
or the double ones, hence the CISD method represents a Full Configuration 
Interaction (FCI) method. In the study of the dimer of ethylene, we expand 
the wave function with a Coupled Cluster Single Double (CCSD) method 
|37I) . This is a numerical tecnique used for describing many-body systems and 
has the property of describing coupled two-body electron correlation effects. 
The CCSD method is an excitation truncated Coupled Cluster (CC) method 
constructed in an exponential approach by including only the desired excitation 
operators (single and double); hence it cannot be considered a FCI- like method. 
However it can describe the interaction of separated molecules better than CISD 
|22| . For both systems, the input is prepared with an Unrestricted Hartree- 
Fock (UHF) calculation. In fact, the UHF description of bond breaking in H2 
gives the proper dissociation products, while the Restricted Hartree-Fock (RHF) 
description of H2 gives unrealistic ones. In the following we show that, in order 
to study the dissociation of H2 molecule, the electron correlation energy must 
be defined by UHF appoximation as |-E? C orr| = \E exac t — Euhf\- In fact, at 
short internuclear distances the RHF and UHF wave functions are identical but 
at large distances, outside the bond region, only UHF reproduces the correct 



10 



progress of correlation energy of H2 that must be zero when the two atoms 
are distant and each electron cannot interfere with the other (Figure 4). The 
correlation energy in Eq. (|13[1 is defined in terms of a complete one-electron 
basis. In practice, however, an incomplete basis must be used for the calculation 
of the correlation energy. The term correlation energy is that used to denote 
the energy obtained from Eq. i|13|) in a given one-electron basis. 

The correlation energy usually increases in magnitude with the size of the 
orbital basis, since a small basis does not have the flexibility required for an 
accurate representation of correlation effects. 

Figure 4 : The figure re- 
ports the electron correlation 
as a function of the inter- 
nuclear distance R of the 
Hi molecule using Gaussian 
basis set 3-21 G with pack- 
age Gaussian. The compar- 
ison between correlation en- 
ergy calculated as \E corr \ — 
\E exact - E RHF \ and calcu- 
lated as \E corr \ — \E exac t — 
Euhf\ is showed. 

In order to confirm this theory, we analyze four different kinds of basis sets 
known in the literature as 3 - 21G, 6 - 31G, 6 - 31G**, 6 - 31++G** and we 
note that the more increase the bases used the bigger becomes the contribute 
to correlation energy. 

The basis sets 3 — 21G and 6 — 31G use Is, 2s and lp atomic orbitals hence 
their contributions to the correlation energy are approximatively the same. The 
basis set 6 — 31G** instead includes the orbital Id, while 6 — 31++G** uses the 
same orbitals but these are more diffuse. Hence, we can see that the energy 
correlation obtained with these two bases is bigger than the one obtained with 
the first basis sets. It is however important to note that at short internuclear 
distances the correlation energy strongly depends on the size of orbital basis 
whereas at large distances, outside the bond region, all four courses become 
approximatively the same. 

After the above preliminary study of basis sets we decide to use the smallest 
basis set (3— 21G). This choice allows us to reduce drastically the computational 
cost, which constitutes an enormus advantage in the study of other systems, in 
particular more complex molecules. 




6 Conclusions 

In this paper, we have analyzed and elaborated the suggestion in that the 
entanglement can be used as an alternative estimation of the electron correlation 
in quantum chemistry calculations. 
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We have, firstly, compared the degree of entanglement in the dissociation 
process of Hi molecule calculated as the von Neumann entropy of the reduced 
density matrix S(pi), as in \TZ\, with the one calculated by the entropy of the 
density matrix of the compound system S(p). The two behaviours obtained in 
this way are shown to be similar, hence, we have choosen to adopt S(p) for 
a qualitative estimation of entanglement in more complex systems. Then, we 
have verified that the electron correlation energy for a dimer of ethylene is well 
reproduced by the entanglement S(p) for different configurations of the system, 
as it has been shown in Figure 2. 

Analyzing interacting molecules, we have introduced the interaction electron 
entanglement in order to calculate the only entanglement due to the interaction 
of the two molecules without their internal correlations. As it has been shown 
in Figure 3, there are configurations that maximize the degree of entanglement. 
This study of the degree of electronic entanglement is a useful resource for 
quantum computers whose input states are constructed in order to be maximally 
entangled. 



Appendix A 

The aim of this appendix is clarifying how one can calculate the reduced density 
matrix as that one we propose in Eq. ©. In order to semplify the calculations, 
we consider two electrons in a two-levels system. Let us introduce the Hilbert 
spaced = \^ 2 {l)®y 2 {l)]®\^ 2 {2)®,9' 2 {2)\ = ^ 4 (l)(g)^ 4 (2), where & and 

5? represent the orbital and the spin space, respectively; in brackets we denote 
one of the two electrons and the apex represents the dimension of the space. An 
orthonormal basis for each space ^f 4 in the occupation number representation 
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be written as |W) = J2b=i w a,£> c a c b|0)i where |0) is the vacuum state, 

the coefficient ui, accordingly to Pauli exclusion principle, satisfies the following 

requests: < UJa ' b _ UJb ' a and and c are the creation and annihilation op- 

erators of single-particle states, respectively whose action on the vacuum state 
is 

4ct|0, = W , Me{M,3,4) J;Kt> 3-jn.tt 

Let us now consider the four sites instead of the two particles: they are single 
particle sites consequently they have no particle or one. Since our aim is having 
some informations about the n\ we proceed making the partial trace of the 
density operator p with respect to n-i \ and n-i J, that is, 
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Pm=Tr n3 p = Y, „ 2T= o,i ( n 2 T,^2 4 |p|n 2 T,«2 4) = 

"2 1= 0, 1 

= (00|p|00) + (01|p|01) + (10|p|10) + (ll|p|ll). 
Let us calculate each matrix element separately. 

The first one is (00|p|00) = Ea l6 =iE«',6'=i w a,6w:^<00|ctcJ|0)(0|c^Ca'|00) 
where a, 6, a' and b' assume the values 1,2,3,4 and they represent the four 
sites of the two levels system that we are studying; moreover, since we make 
the inner product with (00 1 and |00), the only elements of the summatory that 
have non-zero contribution are 1 1 100) 1 thus they correspond to u>i,2 and w 2 ,i : 
|1100): |m T) ® K 4) -»■ wia 
\ni I) ® \ni t) -> w 2 i. 
Hence (00|p|00> = E^ =1 £*, )6 , =1 w„, b <, ji) ,(00|44|Q)(0| C6 ,c a ,|00). Making 
explicit the sum and using the fact that u>a = 0, it becomes: 

(00|p|00) =a; l! 2^ )2 <00|cl4|0)(0|c 2 ci|00)+a;2,i^ )2 <00|4c t 1 |0)(0|c2Ci|00)+ 

+w li2 w* 1 (00|c t 1 c 2 |0)(0|cic 2 |00)+w 2il w* 1 (00|c 2 c t 1 |0)(0|c 1 c 2 |00). (17) 

The action of the creation and distraction operators on the vacuum state pro- 
duces some states such as ±|ni f n\ [) and ±|ni 4 n\ ]) where the double sign 
depends on the order of the operators; using no more the representation refering 
to fermions but the one refering to the sites all these states can be represented 
with ±|1100) so Eq JT7J) becomes: 

(00|pj00> = wi, 2 c^ 2 (00|1100)(1100|00) - w 2 ,iw^ 2 (00|1100)(1100|00)+ 

-o; l! 2w| )1 (00|1100)(1100|00) + ^2,1^2,1 (00|1100) <1100|00) . (18) 

Using the orthogonality, we obtain: (00|p|00) = 4|wi 2 | 2 |ll)(ll|. Let us analyze 
the second matrix element: 

(01|p|01) = El 6= iEt, 6 ' = i^,K', fc '(01|4c!lO>(0|c6'C a n01). Since we make 
the inner product with (01 1 and |01), the only elements of the summatory that 
have non-zero contribution are 1 1 100) thus they correspond to wi j2 and £^2,1' 
|1001): \m t) <8> \n 2 4) -» w M 
j«2 I) O jm T) -> W41 
|1001): |m I) ® \n 2 i) -> cj 24 
K I) <S> |ni |) -> w 42 . 

Hence (01|p|01) = E a ,6=i,2,4 E a ',6'=i,2,4 Ua,V*l,, v (01|ctc+|0) (0|c 6 ,c a , |01). 
Making explicit the sum above and using the site representation as before, this 
matrix element becomes: (01|p|01) = 4|wi.4| 2 |10)(10|. Now, for the third matrix 

element, that is: (10|p|10) = E q,6=i E^=x u*,V»l;v (10|44l )(°l c fr' c »'l 10 ) 

1 In all these calculations, each vector with four components has the first and the second 
entries corresponding with each two sites of ni level and the third and fourth entries corre- 
sponding with n,2 level. Hence the ket 1 1100) , for example, represents two electrons on the 
level ni and no one on the level rt2. 
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we note that the only vectors that have a non zero contribution in the summar- 
tory after the inner product with (10| and 1 10) are: 
|1010): \m T) O |n 2 T) -» W13 

j«2 T) ® |m T) -> w 3 i 
I0110): jm |) ® |n 2 T) -» ^23 

K T) ® jril I) -> ^32- 

Hence, (10|p|10) = 4|cj 2 ,3| 2 |10)(10|. Finaly, the latest matrix element, is, 
analogously, (ll|p|ll) = 4|w 3i4 | 2 |00)(00| 

Thus, in the basis ||00); |01); |10); p ni reads 



/ 4|^ 3 ,4| 2 \ 

4|w 2i3 | 2 

4|w M | 2 

V 4|u lj2 | 2 j 



(19) 



Equation (|19|) represents the density operator for the level n\ and it is important 
to spend some words on the entryes of the matrix: 0^4 represents the excitations 
on the second level so, in the two level system that we consider, W3 4 is associated 
with the probability that n\ is empty; u>i t 2 is associated with the possibility that 
the two electrons are in ni, the coefficients 0^2,3 and u>\^ are associated with the 
probability that n\ can be occupated by one eletron. This fact allows us to in- 
troduce the following notation: p„ 1)0 = 4|u; 3i4 | 2 ; p„ lj2 = 4|wi :2 | 2 ; p ni ,i = 

4|w 2 3I 2 \ 
' ., 12 , where p ni denoted an empty orbital, p„. 2 denotes 

two accupied orbitals and p ni ,i denotes one electron occupied orbital. In order 
to obtain some informations about one of the two electrons, we have to make 
tha partial trace of p ni respect one of the sites of the level n\ (jii ] , n\ J,): 

/ 4(|c 3 . 4 | 2 + K 3 | 2 ) 

\ 4(Jw M | 2 + K 2 | 2 y 

In terms of CISD expansion coefficients, the transition amplitude u>ij can be 

written as |wi, 2 | 2 = — , I w mI 2 = T' 1^3,4 1 2 = ^f- and |w 2i3 | 2 = — , hence pi 

becames: pP 515 = | ' Cl ' +J Cl < 2 ' ) . Eq © is obtained as a 

V |c | 2 + \C2\ J 

generalization of this model to a 771-level system. 
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